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ABSTRACT 

Using a Monte Carlo method, we derive spectra arising from Comptonization 
taking place close to a Kerr black hole. We consider a model consisting of a hot 
thermal corona Comptonizing seed photons emitted by a cold accretion disc. We find 
that general relativistic effects are crucial for the emerging spectra in models which 
involve significant contribution of radiation produced in the black hole ergosphere. 
Due to this contribution, spectra of hard X-ray emission produced in the vicinity of 
a rapidly rotating black hole strongly depend on the inclination of the line of sight, 
with larger inclinations corresponding to harder spectra. Remarkably, such anisotropy 
could be responsible for properties of the X-ray spectra of Seyfert galaxies, which 
appear to be intrinsically harder in type 2 objects than in type 1, as reported recently. 
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1 INTRODUCTION 

A large fraction of the energy liberated by accretion onto 
black holes is released in an optically thin plasma, which 
produces hard X-rays by Compton upscattering of soft pho- 
tons. In some cases this hard X-ray source must be located 
in the immediate vicinity of a black hole, as indicated by Fe 
Kq emission line profiles revealed by Chandra and XMM- 
Newton (e.g. Wilms et al. 2001; Miller et al. 2002). If those 
lines are produced at inner accretion discs, their steep radial 
emissivities require illumination by hard X-rays to be cen- 
trally concentrated within a few gravitational radii, which 
in turn implies that the X-rays are generated close to the 
black hole horizon. In such cases the X-ray spectra should 
be affected by Doppler and gravitational energy shifts as 
well as ray bending effects. The influence of strong gravity 
on Compton scattering expresses itself in an extreme way 
in the ergosphere, where scatterings can extract rotational 
energy of a black hole (Piran & Shaham 1977a). 

Despite the expected importance of these effects, there 
has been little work done to investigate their impact on 
Comptonization spectra emerging from black hole accretion 
flows. Piran & Shaham (1977b) considered production of 7- 
rays in the black hole ergosphere by Compton upscattering 
of X-ray photons. However, their results do not have a form 
suitable to compare them with current observations of black 
hole systems. We find also that their conclusions are par- 
tially incorrect. More recently, Comptonization in advection- 
dominated accretion flows around Kerr black holes has been 
analyzed by Kurpiewski & Jaroszyhski (1999) and Manmoto 



(2000). However, emission from such flows is dominated by 
radiation produced in regions rather distant from a black 
hole, where effects of strong gravity are not significant, as 
we find in this paper. Laurent & Titarchuk (1999) derived 
spectra of Comptonization on a relatively cold plasma mov- 
ing with relativistic velocities close to a black hole. However, 
these bulk motion Comptonization models have assumed so 
far only the Schwarzschild metric, omitting processes tak- 
ing place in the ergosphere, where strong gravity effects are 
most prominent. 

In this paper we study formation of Comptonization 
spectra using a fully general relativistic description of pho- 
ton transfer and Compton scatterings in a plasma located 
close to a Kerr black hole. For our investigation we choose 
a model in which seed photons are emitted by a cold disc, 
while Comptonization takes place in a hot thermal corona. 
We derive the cold disc emission and the Comptonization 
spectra by following individual photon trajectories. We use 
a Monte Carlo method to simulate both the thermal emis- 
sion of the disc and Compton scatterings. 

We take into account two geometries of a disc-corona 
system which are popular in modeling the central part of 
black hole accretion flows. Firstly, we consider a model in- 
volving inner corona replacing a cold disc in the central re- 
gion. This geometry is similar to solutions with a hot opti- 
cally thin flow existing in the inner part of a cold disc (e.g. 
Shapiro, Lightman & Eardley 1976; Narayan & Yi 1995; 
Esin, McClintock & Narayan 1997), although we assume a 
much lower distance of transition between the disc and the 
inner corona than in standard versions of this solution. Sec- 
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ondly, wo consider a corona surrounding a cold disc, corre- 
sponding to models of a hot dissipative corona formed on 
top of the disc surface (e.g. Liang & Price 1977; Haardt & 
Maraschi 1991). 

In Section 2 we describe details of our model. In Section 
3 we analyze relevant relativistic effects. Examples of our 
model spectra are shown in Section 4. We discuss our results 
in Section 5 and summarize our conclusions in Section 6. 



2 MODEL DESCRIPTION 

2.1 Geometry and model parameters 

We consider a Kerr black hole accreting matter at a mass ac- 
cretion rate M . The black hole is characterized by its mass, 
M, and angular momentum, J. We use the Boyer-Lindquist 
(BL) coordinate system x % = (t, R,9, <j>). The metric ten- 
sor components, gij, are given, e.g., in Bardeen, Press & 
Teukolsky (1972). We also make use of locally non-rotating 
(LNR) frames (Bardeen et al. 1972), which are particularly 
convenient for studying processes in the Kerr geometry. 

The following dimensionless parameters are used in the 
paper 

_ R_ - _ ct_ . _ M J o _d0 

where M E dd = 47rGMm p /(<7Tc) is the Eddington accretion 
rate and R g — GM/c 2 is the gravitational radius. Inclination 
of the rotation axis of the black hole to the line of site is given 
by /i bs = cos S b s . We focus on black holes characterized 
by the highest expected value of the spin parameter, a = 
0.998 (Thorne 1974). Such black holes are referred to in this 
paper as maximally rotating. Processes taking place in the 
ergosphere, which is a region contained within the surface 
given by R CIg = Rg[l + (1 — a 2 cos 2 9) 1 ^ 2 ] but outside the 
event horizon, appear particularly important for our study. 

We assume that the accreting material forms a geomet- 
rically thin, optically thick disc and that a hot coronal region 
exists in the neighborhood of the innermost part of the disc. 
We consider two models of a disc-corona system: 

(i) an inner corona replacing a cold disc inside some transi- 
tion radius, r tr ; 

(ii) a cold disc extending down to the black hole horizon with 
a hot corona formed above the surface of the innermost disc; 
the radius of a region covered by the corona in this model 
is denoted as r c . 

We focus our investigation on the impact of strong grav- 
ity effects on formation of spectra rather than on a pre- 
cise modeling of a disc-corona system. Therefore, we make 
a number of simplifying assumptions in our model. 

We follow standard description for a geometrically thin 
disc (e.g. Novikov & Thorne 1973), assuming that the disc 
rotates in the equatorial plane of the Kerr geometry with 
Keplerian angular velocity (Bardeen et al. 1972), 



and that it emits released energy as blackbody radiation 
with the surface brightness distribution described by the 
standard formula given by Page & Thorne (1974). For a 
disc extending down to the black hole [model (ii)] there is 



no emission below the marginally stable orbit, located at 
r ms = 1.23 for a = 0.998. 

In our model of a corona we assume that hot electrons 
have Maxwellian distribution of energy in the corona rest 
frame; that the corona is isothermal; and that the coronal 
material corotates with the disc (except for model B defined 
below) with no transversal or radial motion, u e — and 
u r = 0. 

While corotation of the hot plasma with the disc seems 
natural for a model with corona formed above the disc sur- 
face, velocity field in the inner corona is an uncertain issue. 
However, hydrodynamical studies typically indicate that in 
the extreme Kerr metric the region where radial velocity 
becomes relativistic is very narrow, which supports our ap- 
proximation of the velocity field by a quasi-Keplerian mo- 
tion. 

In particular, solutions of advection-dominated accre- 
tion flows (Popham & Gammie 1998; Manmoto 2000) indi- 
cate that for a rapidly rotating black hole accretion proceeds 
with a rather low radial velocity even close to the event hori- 
zon (while for slowly rotating black holes rapid acceleration 
of the flow occurs at further distances). Although these re- 
sults are not directly applicable to our model as they are 
based on assumption that a quasi-spherical flow is formed far 
from a black hole, they confirm that fast rotation of a black 
hole stabilizes rotationally supported flows. Specifically, for 
a — 0.998 these models typically yield radial velocity in- 
creasing from 0.3c at r = 2 to 0.7c at r = 1.2. Emission of 
accretion flows with such velocities differs from emission of a 
Keplerian flow mainly by increased fraction of photons cap- 
tured by black hole (see Section 3.2). Then, overestimated 
strength of emission from the ergospheric region is a major 
shortcoming due to neglect of radial velocity in modeling 
spectra of such flows. 

Even stronger support for our model comes from the 
MHD simulations, which do not reveal formation of any 
significant sub-Keplerian component. In simulations with 
a = 0.998 a quasi-Keplerian torus extends deep into the er- 
gosphere and the plunging region is extremely narrow (the 
radial velocity profile on Figure 12 in De Villiers, Hawley & 
Krolik 2003 shows velocity smaller than 0.1c at r ms ). 

On the other hand, our description of the accretion 
flow is not appropriate for investigation of spectra pro- 
duced by strongly sub-Keplerian flows, which are postulated 
in models considering Comptonization dominated by trans- 
fer of kinetic energy from bulk motion of such a flow (e.g. 
Chakrabarti & Titarchuk 1995, Laurent & Titarchuk 1999). 
We make a thorough study of spectra arising from such rel- 
ativistic inflows in a separate paper, extending formalism 
developed here by a proper treatment of the propagation of 
photons in a plasma with high radial velocity. 

Velocity of the corona rest frame relative to the LNR 
frame and corresponding Lorentz factor are related to an- 
gular velocity, Q(r, 9), by (Bardeen et al. 1972) 

sin OA . „ r , -,,-2-1-1/2 

V =-=x£T72^-u J ), r ln =[i-v 2 ] , (2) 

where 

A = r 2 — 2r + a 2 , E = r 2 + a 2 cos 2 0, 

A= (r 2 + a 2 ) 2 -a 2 Asm 2 9, uj = 2ar/A. (3) 
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The rest density of electrons in the corona, n, is de- 
scribed by a dimensionless parameter, 



h — nRe-aT, 



(4) 



where at is the Thomson cross-section. 

The following assumptions on distributions of the rest 
density and angular velocity define three specific models of 
the corona. 

Model A: a spherical corona. This model assumes that the 
distribution of rest density is uniform, n — const, in a sphere 
with radius r tl - for model (i) and r c for model (ii), and that 
the angular velocity is constant on spheres, with 



n(r,0) = n K (r). 



(•») 



Model B: a non-rotating spherical corona. This model 
adopts the same assumptions as model A except for angular 
velocity which is set to be equal to the velocity of the LNR 
observer, f2 = u, yielding V = 0. 

Model C: a slab-like corona. This model assumes that an- 
gular velocity distribution is given by 



Q(r,6) = [a + (rsini 



,3/2-1 



(6) 



and that the corona fills only a region within 12° from 
the equatorial plane. The latter assumption is due to the 
fact that angular velocity given by equation © violates the 
causality condition, gtt + 2g tr j,Q + g^Q 2 < 0, at high lat- 
itudes. In this model we assume also that the rest density 
decreases with radius asnar -1 , which gives approximately 
constant vertical optical depth of the corona. 

Most of our model spectra, in both geometries of a disc- 
corona system, are computed for a spherical corona model A. 
Spectra of discs sandwiched by coronae described by models 
B and C are presented in Section 4.2 to illustrate influence of 
assumed velocity field and density distribution on spectral 
properties. 

We do not attempt to model a region where the flow ac- 
celerates inward to relativistic velocities. As discussed above, 
this region should be very narrow and its emission is very un- 
likely to give any noticeable contribution to observed spec- 
trum. Specifically, we neglect the region inside r m i n = 1.13, 
where approximation by a Keplerian motion would yield ar- 
tificially high orbital velocity (V > 0.6, approaching the 
speed of light at the prograde circular photon orbit located 
at r p h = 1.07). We have checked that spectra presented in 
this paper are weakly (in most cases negligibly) sensitive on 
the choice of any specific value of r m i n close to r ms . The most 
significant change of spectrum due to the change of this pa- 
rameter occurs for edge-on spectrum in Fig. Q^a) , which has 
a strong contribution of emission from deep regions of the 
ergosphere (spectra obtained with r m i n = 1.13 and 1.23 are 
compared in that figure). 

2.2 Photons in the Kerr geometry 

We summarize here relevant properties of photon motion in 
the Kerr metric. A photon trajectory is determined by two 
dimensionless constants of motion 



A 



Lc 



where Ei n f is the photon energy at infinity, Q is the Carter's 
constant, L is the component of angular momentum paral- 



lel to the black hole rotation axis. The following equations 
govern the photon trajectory (e.g. Bardeen et al. 1972) 

dt_A(l-u\) dr _ R 1 / 2 d(9 _ 9 1/2 
dC "~ EA ' dC "~ E ' d( ~~ E ' 



dcj> _ A(E - 2r) 2ar 
dC ~~ EA sin 2 6 + EA' 

where 

R = (r 2 + a 2 - \a) 2 - A [(A - a) 2 + 77] 
9 = X] + cos 2 6 (a 2 - A 2 / sin 2 6) , 



(8) 



(9) 



and C, is an affine parameter. 

While photon motion is simply described in the BL co- 
ordinates, local rest frames are more appropriate for descrip- 
tion of physical processes. Applying the transformation be- 
tween the BL coordinate frame and the LNR frame, given 
by Bardeen et al. (1972), we find energy and momentum of 
a photon measured by the LNR observer 

p p n ^( A V /2 E ^ x r E V /2 

£ ln = £ inf (l-aA)^— J , PW-—— e {^) ■ 



pee) 



^(f) 



P(r) 



(10) 



c vea; 

Then, a Lorentz transformation yields the energy and mo- 
mentum in the rest frame of the plasma. In particular, pho- 
ton energy in the rest frame 

Brest = T ln (E ln - . (11) 

Putting in the above equations velocity of a Keplerian 
disc, = Vr = (^k — lo)A/(ZA 1/2 ), we find direct relations 
(given in a similar form by Bardeen & Cunningham 1973) 
between the constants of motion of a photon, A, rj and £"i n f , 
and its energy and directional angles measured in the disc 
rest frame 



A = 



n 



sin <j) cm sin 6> cm + V K 

! K sin cm sin <9 cm + Y}A 1 / 2 A~ 1 



+ uiV K '' 



(12) 



where E^iec is the energy in the disc rest frame, # cm is the 
polar angle with respect to direction perpendicular to the 
disc surface, and 4> e m is the azimuthal angle, in the disc 
plane, with respect to coordinate r direction. 

For motion in a curved space-time we must integrate 
optical depth along a photon trajectory separately for each 
photon. The increase of the optical depth is given by dr = 
cr(E les t)ndl' , where a is the Klein-Nishina cross-section (av- 
eraged over Maxwellian distribution of electron velocities) , n 
is the rest density of electrons and dl' is the length along the 
photon path measured in the plasma rest frame. We find this 
length using the obvious relation, dl' = cdt' , where dt' is the 
time interval measured in the plasma rest frame which can 
be derived for a given trajectory as follows. Transforming 
the BL coordinate differentials dt and d(j> to the LNR frame 
and then making a Lorentz transformation to the plasma 
rest frame we get 

nlt ' -- r ln [(AZ/A) 1/2 dt 



Re 



-V{A/Y<) 1/2 sm8 (d(f>-u>di)] 



(13) 
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Figure 1. Optical depths along radial trajectories of outgoing 
photons emitted at r = 1.2. Solid, dotted and dashed curves cor- 
respond to trajectories in the equatorial plane of a Kerr black hole 
(a = 0.998). Dot-dashed curve is for trajectory with a constant 
polar angle 9 = (close to the polar axis). Photons propagate 
through a thermal plasma with kT e = 80 keV and with the rest 
density h = 0.4 for the dashed curve and h = 0.2 for the remain- 
ing curves. Dotted curve is for a photon with E[ n f = 50 keV, 
the remaining curves for Ei n { = 0.1 keV. Note that the optical 
depth for a photon with i?i n f = 50 keV is reduced by 50 per cent 
between r = 1.2 and r = 2 (compare dotted and solid curves). 



Then, using equation £IJ, we obtain 



.<x(£ reat ) / AE\ 1 / 2 
n 1 i., I I 

. A sin i 



0"T V A J 



1 + V 



2ar sin 8 \ At 



EAV2 J d( 



-V 



SA 1 / 2 d< 



(14) 



where the derivatives dt/d( and d<^/d£ are given by equation 
©. 

Note that in the ergosphere -E rost is typically a few 
times higher than E[ n { due to the gravitational blueshift. 
Therefore, close to a black hole the Klein-Nishina correc- 
tions can decrease significantly the optical depth even for 
photons which in fiat space-time would be still in the Thom- 
son regime. Fig.0sh.ows optical depths integrated along ra- 
dial trajectories (characterized by A = and pg = 0; the 
latter condition implies r\ = —a? cos 2 8) of photons escaping 
from r = 1.2. Photons propagate through a plasma with the 
same temperature and rest density as assumed for models 
presented in Section 4. Note that the optical depth for a 
photon with Ei n { = 50 keV is decreased by around 50 per 
cent in the ergosphere due to the reduction of the cross sec- 
tion. Then, curvature in spectra may occur at a few tens of 
keV (indeed seen in some spectra presented in Section 4). 

Note also that assumptions underlying a spherical 
corona model yield similar optical depths of the corona in 



different directions with respect to the symmetry axis of the 
system. 



2.3 Monte Carlo simulation 

Our procedure for computing spectra is based on the Monte 
Carlo method and involves calculation of many photon tra- 
jectories. A similar approach to computing spectra of black 
hole accretion flows was applied, e.g., by Laor, Netzer & Pi- 
ran (1990), Laurent & Titarchuk (1999) and Kurpiewski & 
Jaroszyriski (1999). The free parameters of our model are rt r 
or r c (depending on the disc-corona model), a, M, m, h and 
electron temperature, T c . Note that M and rh are relevant 
only for the emission of the cold disc. We trace individual 
photon trajectories originating at the disc surface. The pro- 
cedure for generation of disc photons is described in Section 
2.3.1 below. 

Trajectories of photons reaching the hot plasma are fol- 
lowed, through subsequent scatterings, according to algo- 
rithm described in Section 2.3.2. Photons reaching infinity 
are summed over 10 bins in /j, oos (from 0-0.1 to 0.9-1; fi b B 
given for spectra in Section 4 correspond to the value in 
the middle of the bin). Photon transfer equations are solved 
down to r m i n = 1.13. Photons crossing the sphere of this 
radius are treated as captured by the black hole. 

Spectra presented in Section 4 have been obtained in 
Monte Carlo runs involving at least 10 8 photons scattered 
in the corona. 



2.3.1 Thermal disc emission 

The emission point, r cm , on the disc surface is generated ac- 
cording to radial emissivity, F + , given by Page and Thorne 
(1974). Consecutive photons originating from this r em are 
generated until their summed energy exceeds some fixed 
value, Etot, much higher than a characteristic photon en- 
ergy. When Etot is exceeded we generate next r cm . 

For each photon its energy, -Edisc, is generated from 
the blackbody distribution with the local surface temper- 
ature, T sur f = (F + /a) 1 / 4 . The initial direction in the disc 
rest frame is generated with a uniform distribution in e m- 
The polar angle, # om , is generated according to angular dis- 
tribution of local disc emission, with the specific intensity 
7(cos# cm ) ~ 1 + 2.06# em corresponding to limb darkening 
in electron scattering limit (see, e.g., Gierlinski, Maciolek- 
Niedzwiecki & Ebisawa 2001). Constants of motion, 77, A and 
Ei n f, are found from equation [1 1 2 i . Then, we solve r and 8 
equations of the set of equations JSJ and determine whether 
the photon escapes to infinity (contributing to the spectrum 
of the disc emission) , returns to the disc (we neglected here 
increase of the disc temperature due to these returning pho- 
tons), illuminates the corona or gets captured by the black 
hole. 



2.3.2 Comptonization 

If a photon reaches the coronal region, the optical depth is 
generated according to the decreasing exponential distribu- 
tion P(t) — e~ T . Then, we solve equations JSJ and integrate 
the optical depth, given by equation (11411 . until the trajec- 
tory crosses the horizon (see Section 2.3.3), leaves the corona 
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or reaches the chosen optical distance. In the last case a sim- 
ulation of Compton scattering is performed as follows: 

(i) Photon momentum and energy in the plasma rest frame 
are found, as described in Section 2.2. 

(ii) Electron velocity is generated and scattered photon di- 
rection and energy are derived according to the procedure 
described in Gorecki & Wilczewski (1984). 

(iii) The four-momentum of the scattered photon is Lorentz- 
transformed back to the LNR frame, then equations 1101 are 
inverted for new constants of motion (i?i n f, t], A). 

Photons leaving the corona are treated similarly to 
those emitted from the disc surface - considering their mo- 
tion in the r6 plane we determine whether they escape from 
the system and contribute to the Comptonization spectrum. 



3 COMPTON SCATTERING IN THE KERR 
METRIC 

3.1 Single scattering spectra 

Photons scattered in a hot corona gain energy due to both 
thermal motion of electrons and bulk motion of the corona. 
Furthermore, dynamics of individual scatterings is influ- 
enced by gravitational energy shifts. In this section we de- 
termine the range of distances from a black hole in which 
the gravitational effects are important for the formation 
of Comptonization spectra. Then, we examine qualitatively 
how these effects can affect the emerging spectra. For this 
investigation we derive spectra of a single scattering of mo- 
noenergetic photons, with initial energy E® n{ , on electrons in 
the Kerr geometry. We assume various distributions of elec- 
tron energies to compare effects taking part in formation of 
the spectra. Simulation of Compton scatterings is performed 
according to procedure described in Section 2.3.2. 

Spectra presented in this section are normalized by the 
number of impinging photons, n p h(-Emf) = N(E- m f)/No, 
where No is the number of incident photons and N(Ei a f) 
is the number of photons observed, after scattering, at in- 
finity with energy E ln f. Note that photons captured by a 
black hole do not contribute to the spectrum, therefore 

J n ph {E ini )dE int < 1. 

To illustrate effects of general relativity, separated from 
kinematic and thermal effects, we examine first scattering 
on electrons at rest with respect to the LNR frame. Fig. 
|2I a) shows spectra for scattering of photons, with Ef nf — 30 
keV, on electrons located in the equatorial plane of a Kerr 
black hole. We assume that photons illuminate electrons 
isotropically (as seen in the LNR frame) excluding directions 
which correspond to negative energy, Ef nt < 0, of imping- 
ing photons. Also shown is the spectrum, integrated over 
all scattering angles, for scattering of photons with ener- 
gies Eo = 30 keV on electrons at rest in flat space-time. 
For Eo <C m c c 2 Compton scattering is nearly elastic and 
the spectrum emerging in flat space-time is very narrow. 
On the other hand, scattering in the Kerr space-time pro- 
duces spectra extending remarkably further into both low 
and high energies. The low energy part of such spectra is a 
generic feature of scattering close to a black hole, resulting 
from the gravitational blueshift of incoming photons. Due 
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Figure 2. Angle-averaged spectra for a single scattering of pho- 
tons, with = 30 keV, impinging on electrons at rest (with 
respect to the LNR frame) in the equatorial plane of the maximal 
Kerr geometry, (a) Photons illuminate electrons isotropically (as 
seen in the LNR frame). Scatterings take place at r = 1.4 (dashed 
curve), r = 3 (dot-dashed curve) and r = 6 (dotted curve). Solid 
curve gives the spectrum for a single scattering on electrons at 
rest in flat space-time, integrated over all scattering angles. The 
solid and dotted curves are rescaled by a factor of 0.1. Photons 
emerging from r = 1.4 with energies exceeding 135 keV (indicated 
by the solid vertical line) come from scatterings which extract the 
black hole rotational energy, (b) Dependence of the spectrum for 
scattering at r = 1.4 on the direction of impinging photons. The 
initial direction of photons is given by the following angles (de- 
fined in the text; fi\ n = cos6*i n ): fi\ n = 0, <j>i n = it for dot-dashed 
curve; fn n = 0, <ft n = tt/2 for dotted curve; fn n = 1 (which 
uniquely determines the direction) for solid curve; and fi\ n = 0.2, 
<frl n = 1.917T for dashed curve. Scatterings extracting the rota- 
tional energy take place only for photons impinging from the last 
direction. 



to this blueshift, close to a black hole photons lose in a scat- 
tering, on average, a larger amount of energy than in flat 
space-time. 

More importantly, in the Kerr metric the relation be- 
tween photon energies in the LNR frame and at infinity in- 
volves photon angular momentum, see equation 1101 . There- 
fore, the photon energy after scattering may be higher than 
its energy before scattering, provided that a suitable change 
of angular momentum occurs. Such effect is not allowed in 
either flat space-time or the Schwarzschild metric. The in- 
crease of photon energy can be particularly significant in 
the black hole ergosphere, where Penrose Inverse Compton 
scatterings (Piran & Shaham 1977a), for which a scattered 
photon has higher energy than the sum of the photon and 
electron energies before scattering, are permitted. Photons 
produced in such events form the high energy part of the 
spectrum, as indicated on Fig. |2Ia) for the spectrum from 
r = 1.4, which extends to energies exceeding the initial en- 
ergy by over an order of magnitude. 
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Figure 3. Angle-averaged spectra for a single scattering, at 
r = 1.4, on electrons at rest in the equatorial plane of the 
maximal Kerr geometry. Solid, dotted and dashed curve is for 
E9 f = 30, 300 and 3000 keV, respectively. Photons illuminate 
electrons isotropically, as seen in the LNR frame. 



The necessary condition for a Penrose scattering to take 
place is a very large blueshift of an impinging photon in 
the LNR frame. The blueshift is determined by the angle, 
"^ln, between the photon momentum and the <f) direction in 
the LNR frame. Namely, using equation 1101 . with — 
(Ein/c) cos $i n , we get 



E ]n = £ illf (E4) 1/2 (EA 1/2 + 2ar sin (9 cos $ b 



(15) 



Increasingly high blueshifts correspond to narrower ranges 
of initial directions. E.g., for E ln f = 30 keV and r = 1.4, Ei n 
exceeds 100 keV for 0.45tt < $ hl < 0.59tt, while E ln > 1000 
keV for 0.58-7T < $i n < 0.597T. Note also that strongly 
blueshifted photons have very large negative angular mo- 
mentum (A <C 0) and can emerge only from a radiative 
process taking place in the ergosphere. 

Clearly, Penrose scatterings are permitted only for cer- 
tain initial directions of a photon. Fig.[5f b) illustrates depen- 
dence of the scattering spectrum on initial direction. Note 
that while the blueshift is determined by $i n , different initial 
directions with the same <J>i n give rise to different scatter- 
ing spectra at infinity. The spectra in Fig. |5fb) correspond 
to initial directions determined in the LNR frame by the 
polar angle, #i n , between the photon initial direction and 
coordinate ^-direction, and the azimuthal angle, cj}\ n , in the 
r<j> plane, with respect to the r-direction. These angles are 
related to $i n by cos$i n = sin <^i n sin #i n . 

Considering possibility of generating high energy pho- 
tons due to Penrose processes in black hole accretion flows, 
we note that a fine tuning of the initial direction is required 
for the most efficient scatterings and therefore probability 
of such events is low. Moreover, even if a photon propagates 
with sufficiently high blueshift in a local rest frame, proba- 
bility of a scattering is strongly reduced due to the decline 
of the Klein-Nishina cross-section. 

Fig.|3]shows dependence of the scattering spectra on ini- 
tial energy of photons. The fractional loss of photon energy, 
in the LNR frame, increases with increasing Ef nt . There- 
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Figure 4. A spectrum for a single scattering of photons, with 
Eo = 30 keV, on electrons with a Maxwellian distribution of ve- 
locities, kT c = 80 keV, in flat space-time (given on both figures 
by the solid curves) compared with angle-averaged spectra for a 
single scattering of photons, with E9 . = 30 keV, on electrons 
in the equatorial plane of the Kerr geometry, (a) Scattering at 
r = 1.4. Dashed curve is for scattering on cold electrons at rest 
in the LNR frame. Dot-dashed curve is for scattering on thermal 
electrons, with kT c = 80 keV, at rest in the LNR frame. Dotted 
curve is for scattering on thermal electrons (kT c = 80 keV) on a 
Kcplcrian orbit, (b) Scattering at r = 6 for a = (dashed curve) 
and a = 0.998 (dotted curve). Target electrons have thermal en- 
ergies (kT c = 80 keV) and rotate on a Keplerian orbit. Note that 
difference between spectra emerging in the Schwarzschild and ex- 
treme Kerr geometries is negligible. 



fore, for higher energy photons the gravitational effects are 
relatively less efficient in increasing photon energies. 

Outside the ergosphere, the strong gravity effects are 
not so noticeable, as illustrated on Fig. Ufa) . Photons scat- 
tered at r = 3 achieve energies exceeding the initial energy 
by up to 50 per cent and deviation of spectrum emerging 
from scatterings at r — 6 from the flat space-time scat- 
tering spectrum is insignificant. Therefore, outside the er- 
gosphere other effects should be dominant in formation of 
Comptonization spectra. 

Fig. 0| compares efficiency of raising the energy due 
to gravitational shifts with boosting of photon energies re- 
sulting from transfer of thermal and bulk motion energies 
of electrons. Note that the probability that a photon with 
E? n! = 30 keV emerges with energy exceeding 100 keV after 
scattering on a cold electron at r = 1.4 is equal to that corre- 
sponding to a single scattering in a semi-relativistic thermal 
plasma in flat space-time. Velocity of the circular motion re- 
mains relatively low even in the deep ergosphere (ri n < 1.3). 
However, combined effects of strong gravity, thermal and 
bulk motions may result in a large gain of energy in a single 
scattering, as illustrated by the dotted curve on Fig. Ufa). 
Then, an optically thin plasma in the ergosphere is able to 
generate photons with energies a few times higher than ther- 
mal energy of electrons. 

Outside the ergosphere, thermal effects completely 
dominate. Spectrum emerging from thermal electrons in cir- 
cular orbit at r — 6 differs only slightly from spectrum pro- 
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Figure 5. Spectra for a single scattering of monoenergetic pho- 
tons, with E? n{ = 30 keV, illuminating isotropically (in the LNR 
frame) electrons at r = 1.4 in the extreme Kerr metric. Dashed, 
solid and dotted curves give spectra of photons scattered into 
^obs = — 0.2, 0.4 — 0.6 and 0.8 — 1, respectively, (a) Spectra 
for cold electrons at rest in the equatorial plane; angle averaged 
spectrum is given by the dashed curve on Fig. 0a). (b) Spectra 
for electrons, on a Keplerian equatorial orbit, with a Maxwellian 
distribution of velocities (kT c = 80 keV) in the rest frame; angle 
averaged spectrum is given by the dotted curve on Fig. Q[a). (c) 
Spectra for electrons with a Maxwellian distribution of velocities 
(fcT e = 80 keV) in the rest frame. The rest of electrons rotates on 
a circular orbit, at 9 = 7r/3, with angular velocity given by equa- 
tion JFJ- (d) Same as in Fig. 0c) but the spectra are formed by 
photons which do not cross the equatorial plane after scattering. 
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Figure 6. Angle-averaged spectra for a single scattering of mo- 
noenergetic photons, with E? ni = 30 keV, on electrons with var- 
ious velocities: (v r ,v'f') = (0,0.56c) for thinner solid curve [Ke- 
plerian flow, same as dotted curve on Fig. 4(a)]; (—0.3c, 0.5c) 
for dashed curve; (—0.7c, 0.36c) for dot-dashed curve; (— 0.9c, 0) 
for thicker solid curve; and (—0.7c, 0.56c) for dotted curve. Scat- 
terings take place at r = 1.4 in the equatorial plane. Isotropic 
illumination (in the LNR frame) of electrons and Maxwellian dis- 
tribution (kT e = 80 keV) of electron velocities in the rest frame 
are assumed. 



models with an untruncated disc such photons are depleted 
from emerging spectra as they interact with the disc matter. 

We find only a very weak anisotropy of spectra emerging 
from r = 6, regardless of the speed of black hole rotation. 
Then, we conclude that spectra of Comptonization in the 
Schwarzschild metric as well as those emerging at distances 
exceeding a few R g in the Kerr metric should be correctly 
approximated by flat space-time Comptonization codes. 



duced by the same thermal population of electrons at rest 
in flat space-time [see Fig. Bib)]. Moreover, the difference 
between the spectra produced, at r = 6, in the vicinity of a 
maximally rotating and a non-rotating black hole is negligi- 
ble. 

Fig. illustrates dependence of spectra emerging from 
Compton scattering in the ergosphere on the inclination an- 
gle of a distant observer. Most photons emerging from the 
equatorial plane are observed from directions close to that 
plane [see Fig. EI a )] due to the bending of photon trajec- 
tories, which is the well known effect of the Kerr geometry. 
Moreover, this gravitational focusing affects primarily high 
energy photons. As a result, spectra observed edge-on are 
much harder than those observed face-on. Isotropic distri- 
bution of thermal velocities of electrons does not diminish 
this anisotropy significantly, on the other hand circular mo- 
tion further collimates radiation along the equatorial plane 
[see Fig.EJb)]. 

Slightly more isotropic spectra come from Compton 
scattering taking place off the equatorial plane [see Fig. 
Hfc)]. Moreover, trajectories of a significant part of photons 
cross the equatorial plane, as illustrated by Fig. |3d). In 



3.2 Radial velocity 

Fig. |S] shows spectra for a single scattering, obtained sim- 
ilarly to spectra shown in Section 3.1, but for various ve- 
locity fields of Comptonizing electrons. The velocity field is 
parametrized by two components of velocity measured in the 
LNR frame: the orbital velocity, %r, and the radial velocity, 
v r (= A _1 A 1 / 2 « r /u*; Bardeen et al. 1972). 

The dashed and dot-dashed curves show spectra for val- 
ues of and v r typical for solutions of advection-dominated 
accretion in the ergosphere (see, e.g., Figure 1 in Popham & 
Gammie 1998). Clearly, scattering in a Keplerian flow (the 
thinner solid curve) results in the highest fraction of escap- 
ing photons. On the other hand, highest energies achieved 
by scattered photons are similar in both models. 

We emphasize that relatively high (Keplerian) orbital 
velocity is essential for spectra emerging from the innermost 
region of accretion flow, as it collimates emission into direc- 
tions for which escape to a distant observer is possible. The 
thicker solid and dot-dashed curves illustrate how decrease 
of orbital velocity affects the spectrum. Although the de- 
crease of is connected with increase of v r , yielding more 
efficient Doppler boosting (Tin = 1.6 and 2.3 for the dot- 
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dashed and solid curve, respectively), it reduces efficiency of 
generating high energy photons able to escape. 

Spectrum for the only velocity field producing photons 
more energetic than a Keplerian flow (due to both a high 
Lorentz factor, T\ n = 2.3, and a Keplerian orbital velocity) 
is shown by the dotted curve. We emphasize, however, that 
such a velocity field is very unlikely to describe a real flow 
as total specific energy in such a flow would exceed by 30 
per cent its rest energy, while the opposite is expected in a 
dissipative flow. 

We conclude that relativistic radial velocities, which can 
significantly increase the Doppler boosting, are allowed only 
in sub-Keplerian accretion flows, emission of which, however, 
is beamed toward the event horizon and lost. 



3.3 High energy cut-off 

In flat space-time, repeated scatterings of soft photons in 
a mildly-relativistic thermal plasma lead to formation of a 
power- law spectrum extending up to kT c . Photons scattered 
to higher energies achieve Wien equilibrium with electrons 
and the spectrum has a roughly exponential turnover re- 
flecting the thermal distribution of electron velocities. Close 
to a black hole this property is disturbed by energy shifts 
between consecutive scatterings. In particular, for photons 
blueshifted to energies significantly exceeding kT c , thermal 
motions of electrons become less significant and the spec- 
trum is formed primarily due to the recoil effect. This effect 
distorts the high energy cut-off from the shape characteristic 
for thermal spectra, as we illustrate in this section. 

Due to a non-local nature of Comptonization, analysis 
of formation of Comptonization spectra in a curved space- 
time is a very complex issue. Below we apply a simpli- 
fied treatment to derive spectra of consecutive scatterings 
of soft seed photons off thermal electrons, with kT c — 80 
keV, which are in circular motion corresponding to model 
A. We assume that all scatterings take place in a fixed point 
(r — 1.4 and 9 = 0.457r) and that in each scattering pho- 
tons impinge from the same direction, given by (j>i n = 7r and 
jttin = 0. For each photon its initial energy, Ei n f, is generated 
from a blackbody distribution with fcTbb = 100 eV. Then, 
the following steps are repeated to simulate four subsequent 
scatterings of this photon: 

(i) rest energy is found [equations ffllOB and llH l - the as- 
sumed initial direction yields E last = 6£7i n f ; 

(ii) simulation of Compton scattering in the rest frame is 
performed; 

(iii) energy of scattered photon, i?i n f, is found and equations 
of motion are solved; photons escaping to infinity contribute 
to spectrum of given scattering order. 

Figs. [7] and [8] show spectra at infinity and in the rest frame, 
respectively, for second and fourth scattering orders. 

The above procedure roughly approximates effects of 
photon transfer in the Kerr metric. Step (i) accounts for 
significant blueshifts experienced by photons in the immedi- 
ate vicinity of the event horizon. In turn, the assumed initial 
direction accounts for the fact that bending of photon tra- 
jectories results in very anisotropic radiation field close to a 
black hole. In simulations including full treatment of photon 
propagation, presented in Section 4, we find that initial di- 
rections of photons scattered close to the event horizon are 
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Figure 7. Solid curves show spectra emerging from second (a) 
and fourth (b) scattering order, observed by a distant observer 
at /^obs = 0.55. Scatterings take place at r = 1.4 and 8 = 0.45-7T 
(see text for details). Dotted curves show spectra, as a function 
of photon energy at infinity, before scattering. All the spectra are 
normalized to unity in maximum. 
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Figure 8. Figure illustrates increasing role of the recoil effect 
in formation of spectra emerging from higher scattering orders. 
Solid curves show angle-averaged spectra, in the rest frame, of 
photons emerging in second (a) and fourth (b) scattering order. 
Dashed curves show spectra of photons before scattering. 

concentrated along the inward radial direction (given by the 
assumed values of (j>i n and /Ui n )- 

This anisotropic illumination is not important in for- 
mation of spectra emerging from lower scattering orders, for 
which the change of photon energies in the rest frame is 
dominated by energy gain from Doppler effect. In such case, 
photons scattered in various directions have similar spectra 
and the relativistic effects result in only a slight redshift and 
smearing of the observed spectrum. On the other hand, for 
higher scattering orders the angular distribution of illumi- 
nating photons becomes important, as spectra formed due 
to the recoil effect depend on the scattering angle and then 
photons scattered into various directions are subjected to 
different transfer effects, as described below. 

Due to high initial blueshift and moderate redshifts of 
photons escaping after scattering (typical redshift of escap- 
ing photons is 0.5 < g < 1), in lower scattering orders 
photon energies, E^, amplify by a large factor with each 
scattering [see Fig. |7fa)]. even though in the rest frame the 
fractional energy transfer is small [Fig. |SJa)]. 

In higher scattering orders photons lose energy, in the 
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rest frame, forming spectrum shown by the solid curve on 
Fig. Hfb) . Spectrum of photons observed by a distant ob- 
server [solid curve on Fig. E£b)] is also affected by effects of 
photon transfer in the Kerr space-time. In particular, scat- 
tered photons with -E ros t > 1 MeV emerge from forward scat- 
terings and they weakly contribute to the observed spectrum 
as most of them are captured by the black hole. Photons ob- 
served with _Bi n f < 200 keV come from backward scatterings 
and they are additionally redshifted by g < 0.5. 

The recoil and transfer effects give rise to spectrum dis- 
torted from the typical thermal shape, specifically, a flat- 
tened part occurs between 200 and 400 keV in the spec- 
trum on Fig. Efb) . Contribution of higher scattering orders 
to Comptonization spectra should result in similar distor- 
tion of the high energy cut-off, although decreased scattering 
probability resulting from the (necessarily high) blueshift 
diminishes contribution of photons scattered in this regime. 
Spectral breaks due to this effect are indeed seen in spectra 
on Fig. G3» above f 00 keV. 



4 THERMAL COMPTONIZATION SPECTRA 

Emergent spectra are shown in Figs. I5IITUI and 1121151 Each 
figure presents spectra for three inclinations, ^i b s = 0.85, 
0.55 and 0.15. The spectra have been obtained for stellar 
mass, M = IOMq, and supermassive, M = W 7 Mq, black 
holes. For a stellar black hole we assume a distance of d = 5 
kpc and for a supermassive one d = 5 Mpc. All models 
assume a maximally rotating black hole and rh = 0.5. Un- 
less specified otherwise, we assume that corona is spheri- 
cal (model A) and that the electron temperature kT a = 80 
keV. Motivated by results of the previous section we focus 
on models with coronae located within a few gravitational 
radii. 
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Figure 9. Spectra emerging from a hot inner corona surrounded 
by a cold accretion disc. The model assumes a maximally rotating 
black hole with a mass M = 10 7 Mq, at a distance of d = 5 Mpc, 
accreting at rh = 0.5. A uniform, spherical (model A, see Section 
2.1) and isothermal (kT c = 80 keV) corona is assumed to replace 
the disc within rt r = 2 (a) and 4 (b). The rest density n = 0.4 
and 0.2 for rtr = 2 and 4, respectively. The spectrum of thermal 
emission from the cold disc and the Comptonized component are 
shown separately. Solid, dashed and dotted curves correspond to 
^obs = 0.85, 0.55 and 0.15, respectively. The dot-dashed curve 
in Fig. 9(a) shows contribution, to spectrum observed at /^ b s = 
0.85, of photons which experienced scattering inside r = 1.4. The 
thinner dotted curve in Fig. 9(a) shows spectrum for r m i n = 1.23 
(see Section 2.1). 



4.1 Inner corona 

Figs. |U] and I1UI show spectra of models with a hot inner 
corona replacing the disc within rtr = 2 and 4. We assume 
values of the rest density yielding a similar Thomson opti- 
cal depth of the corona in both cases (see Fig. 0, namely 
h — 0.4 for r tr = 2 and h = 0.2 for r tr = 4. Two main spec- 
tral components, emitted by the accretion flow, are shown 
separately. Thermal emission of the outer disc, for the as- 
sumed rh, has a peak around 0.1 keV and at a few keV 
for M — 10 7 Mq and IOMq, respectively. Comptonization 
spectra emerging in models with r tr — 2 and 4 differ signif- 
icantly, although the hot plasma parameters are similar. 

For = 2 virtually bulk of Comptonized radiation 
comes from the ergosphere. The luminosity, slope and cut- 
off energy of the spectrum strongly depend on the inclination 
of the line of sight. This is an obvious consequence of bend- 
ing to the equatorial plane, affecting trajectories of photons 
emerging from the ergosphere (see Fig. . 

For photons observed close to the disc plane amplifica- 
tion of energy in the first scattering is so large that their 
energies significantly exceed the energy of disc photons. As 
a result, a low energy cut-off of a Comptonization spectrum 
may be clearly visible in the total spectrum [e.g., note deficit 
of photons below 1 keV in the spectrum for /i obs = 0.15 on 
Fig. Ha)]. 
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Figure 10. Same as in Fig. |5]but with M = IOMq and d = 5 
kpc. 



The supermassive and stellar black hole models produce 
slightly different spectra in the hard X-ray range. This prop- 
erty is not obvious as in our model a change of M affects 
directly only the temperature of the outer disk and hence the 
seed photons energy. The dependence of hard X-ray spectra 
on black hole mass is due to the fact that photon trans- 
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Figure 12. Dependence of Comptonization spectra on the hot 
plasma parameters in a model with rtr = 2 and M = XOMq. The 
model parameters, except for h and kT e , have the same values 
as in Fig, 1101 a), (a) h = 0.2, kT c = 80 keV; for a comparison, 
Comptonization spectrum with h = 0.4 [same as in Fig. HUf all. 
for /i bs = 0.55, is shown by the thinner dashed curve (b) h = 0.4, 
kT e = 10 keV. 



E (keV) 

Figure 11. Solid curves show the Comptonization spectra of a 
model with M = 10 7 Mq and rt r = 4 [same as in Fig. I^Jb)] for 
A'obs = 0.15 (a) and 0.85 (b). Dotted and dashed curves show con- 
tributions of photons which experienced at least one scattering in 
the ergosphere and photons scattered only outside the ergosphere, 
respectively. Dot-dashed curves show contribution of photons for 
which the last scattering took place in the ergosphere. 

fer in the corona is affected by gravitational energy shifts 
(see Section 2.2). As a result, the corona is optically thin 
(t<0.5) for photons with energies of a few tens of keV, 
emerging from the first scattering order in the stellar black 
hole model, but it is moderately thick (r > 1) for photons be- 
low 10 keV, which arise from first few orders of scattering in 
the supermassive black hole model. Therefore, higher energy 
density of the radiation field establishes in the innermost 
region (where photons achieve highest energies due to grav- 
itational effects) of the corona surrounding a supermassive 
black hole. Contribution of high energy photons generated 
in this region results in the shape of the high energy cut-off 
corresponding to the predominant role of the recoil effect 
(Section 3.3) in the model with M = W 7 Mq. Also, stronger 
contribution of photons upscattered in the innermost region 
[shown by dot-dashed curves on Figs. |5Ja) and HUf all. but 
then scattered in the outer region of the ergosphere into di- 
rections corresponding to low inclination, results in harder 
spectra (above a few tens of keV) observed at /i b s = 0.85 
in the model with M = 1O 7 M0. 

For r tI = 4 most scatterings take place in a region 
where standard thermal Comptonization effects dominate. 
However, the slope and luminosity of the Comptonization 
component remain strongly dependent on inclination. This 
anisotropy is entirely caused by different magnitude of the 
contribution of radiation produced in the ergosphere to spec- 
tra observed from various directions (see Fig. IllH . Photons 
which do not undergo scatterings in the ergosphere form 
spectra with the photon spectral index F(l — 80 keV) ~ 2.1. 
A spectrum with such F would result from thermal Comp- 



tonization, with kT e = 80 keV and r = 0.8(= hrtx), in 
flat space-time. This supports our conclusion that the only 
remarkable influence of the space-time metric on the Comp- 
tonization spectrum occurs for radiation emerging from the 
ergosphere. 

Note also that ergospheric emission is significantly 
isotropized due to scattering in the outer regions of the 
corona. Therefore, contribution of radiation from the ergo- 
sphere to spectrum observed at /i b s = 0.85 is significant in 
model with rtr = 4, although photons escaping directly from 
the ergosphere contribute negligibly to this spectrum. 

As noted in Section 2.2, a spectral curvature may oc- 
cur at a few tens of keV due to the decline of the Klein- 
Nishina cross-section. Indeed, a spectral break between 10 
and 20 keV, related to this effect, is seen on Fig. HHh ). Fur- 
thermore, emission from the ergosphere has a significantly 
harder spectrum than emission from the outer regions and 
their superposition may result in an upturn of the spectrum, 
such as seen on Fig.lTTTbl around 1 keV. Then, a power-law 
approximation of the Comptonization spectrum may be not 
adequate for some regions of the parameter space. 

Our model does not involve an explicit description of 
the corona heating mechanism. However, the energy trans- 
ferred to photons, measured in the corona rest frame, deter- 
mines spatial distribution of heating needed to maintain the 
assumed uniform distribution of temperature. In particular, 
models with rtr = 4 require that around 40 per cent of the 
heating power should be provided to the plasma within the 
ergosphere. Roughly 30 per cent of this energy is lost due to 
capturing of photons by the black hole. However, photons 
upscattered within the ergosphere give a major contribution 
to the total luminosity of Comptonized radiation, ranging 
from 50 per cent at low inclinations to over 70 per cent in 
radiation observed edge-on. 

Hardening of the Comptonization component, with in- 
creasing inclination, is accompanied by weakening of the disc 
component. Thermal emission from the innermost region of 
the disc, within several R s , is gravitationally focused to di- 
rections close to the equatorial plane (Cunningham 1975; 
due to this effect, emission from the disc truncated at rtr = 2 
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Figure 13. Spectra emerging from a spherical corona (model 
A) surrounding the innermost disk surface with (a) r c = 2 and 
n = 0.4; (b) r c = 4 and n = 0.2. Other parameters have values 
indicated in Fig. The dot-dashed curve on Fig. USThl shows 
contribution of photons scattered in the ergosphere to spectrum 
observed at /i b s = 0.15. 
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Figure 14. Same as in Fig. [13] but with M = IOMq and d = 5 
kpc. 



is much stronger, at high inclinations, than emission from 
the disc with r tr = 4, see Figs. 151 and 1101 . However, angu- 
lar distribution of emission from more distant regions of the 
disc is mostly due to the reduction of the projected area, 
with the observed flux given by -Fobs oc /i b s at r > 1. Then, 
the ratio of the Comptonization component to the disc emis- 
sion luminosities changes strongly with inclination. E.g., for 
rtr = 2 the Comptonization component contributes 65 per 
cent to the total luminosity at /j, b s = 0.15 but only 5 per 
cent at /z bs = 0.85. 

Fig. ^| illustrates dependence of the spectrum of ergo- 
spheric emission on electron temperature and rest density in 
the corona. Remarkably, a relatively cold plasma (kT c = 10 
keV) is able to produce spectra extending beyond 100 keV, 
especially at high inclinations [see Fig. I12f b'll . 
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Figure 15. (a) Same as in Fig. 1141 b~> (r c = 4) but for a non- 
rotating spherical corona (model B). (b) Similar as in Fig. 1141 a) 
(r c = 2) but for a corona described by model C with h = 4/r. 



4.2 Corona above the disc surface 

Figs. 1131 and 1141 show spectra of models involving an opti- 
cally thick disc extending down to the event horizon, with a 
spherical corona surrounding the disc surface within r c = 2 
and 4. The parameters of the corona are similar to these 
assumed above for the inner corona models, namely n = 0.4 
for r c = 2 and n = 0.2 for r c = 4; kT c = 80 keV. The in- 
nermost disc immersed in the corona provides a very strong 
supply of seed photons. However, the Comptonization spec- 
tra are softer than spectra emerging from inner coronae with 
the same parameters, which is due to the fact that majority 
of photons upscattered in the ergosphere is absorbed by the 
disc [see Fig. OJc) (d)] . 

In models with r c = 4 the Comptonization spectra are 
dominated by emission from regions outside the ergosphere 
and similar slopes characterize spectra observed from vari- 
ous directions. The ergospheric component [shown on Figs. 
I13f b) and I14f b~)l contributes up to 25 per cent to the lu- 
minosity of Comptonization radiation. However, luminosity 
and cut-off energies of the spectra increase for higher incli- 
nations, which is due to the circular motion of the corona. 

To confirm this last conclusion we show on Fig. I15f al 
spectrum of a non-rotating corona (model B). The Comp- 
tonization spectra are again dominated by emission gener- 
ated outside the ergosphere, which is weakly affected by the 
ray bending or gravitational shifts of energies. As a result, 
anisotropy of the Comptonization component is negligible. 

Finally, we consider sensitiveness of the properties of er- 
gospheric Comptonization spectra on angular velocity distri- 
bution in a corona. Fig. I15f b'l shows spectra emerging from 
a corona described by model C with n — 4/r. The corona 
fills a region within 12° of the disc plane and covers the 
disc surface within r c = 2. In this model all scatterings take 
place close to the equatorial plane in a plasma rotating with 
V ~ 0.6, which results in slightly different angular distribu- 
tion of Comptonized photons with respect to the spherical 
corona model. In the latter model a large number scatter- 
ings take place in a slowly rotating plasma at high latitudes. 
However, the basic properties of the spectra appear not af- 
fected by the change of velocity field and geometrical shape 
of the corona. 
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5 DISCUSSION 

5.1 Comparison with similar studies 

The first study of the impact of the Kerr metric on spec- 
tra generated by a black hole accretion flow was performed 
by Cunningham (1975). He pointed out that most of the 
thermal emission from innermost regions of a disc lying in 
the equatorial plane is gravitationally focused into directions 
close to the disc plane. We find that a similar focusing affects 
also radiation arising from Comptonization in the vicinity of 
a Kerr black hole. 

Piran & Shaham (1977a,b) investigated Compton scat- 
tering in the Kerr metric. They pointed out Compton scat- 
tering, taking place in the ergosphere, as the most astrophys- 
ically feasible Penrose process. As an application of such pro- 
cesses they considered production of 7-rays by a hot plasma 
in the ergosphere. 

As we discuss in Section 3, the gravitational effects may 
indeed significantly increase the energy of a photon scattered 
in the ergosphere. However, probability of the most efficient 
processes is extremely low and they appear not important 
for formation of spectra. Piran & Shaham (1977b) point out 
dependence of the emerging spectra on observer's inclina- 
tion, in agreement with our results. On the other hand, they 
concluded that thermal Comptonization in the ergosphere 
gives rise to spectra with a power-law shape extending into 
MeV range without any exponential-like cut-off. We do not 
confirm this last conclusion as our spectra, derived for simi- 
lar electron temperatures, have generic thermal-like cut-offs. 
We suspect that this discrepancy results from a rather poor 
statistical quality of spectra presented by Piran & Shaham 
(1977b). Their spectra were obtained with a typical number 
of 10 3 photons. On the other hand, we find that at least 
10 6 photons scattered in a corona is needed in a simulation 
to correctly derive the high energy part of angle-averaged 
spectrum. 

Thermal Comptonization in two-temperature 
advection-dominated flows in the Kerr metric was studied 
by Kurpiewski & Jaroszyhski (1999) and Manmoto (2000). 
However, their models assumed a hot flow extending out 
to ~ 1000-R 9 . Therefore, their spectra arise mostly due to 
thermal processes taking place far from a black hole, where 
gravitational field is not important for radiative processes. 
Then, these authors did not notice any direct impact 
of the gravitational field metric on a Compton effect. A 
weak dependence of their spectra on the black hole spin 
parameter is due to changes in velocity field, density and 
electron temperature profiles of the accretion flow. 



5.2 Applicability to black hole systems 

Clearly, effects considered in this paper are important for 
the emerging spectra only if a substantial part of hard X- 
rays comes from the ergosphere. On theoretical ground, sce- 
nario with strong emission from the ergosphere, in case of 
accretion onto a rapidly rotating black hole, is supported by 
preliminary results of general relativistic MHD simulations 
(see below) as well as by the standard theory of Keplerian a- 
discs. The latter predicts that 17 per cent of the gravitational 
power is dissipated within the ergosphere (for a = 0.998). 
This implies that significant emission of hard X-rays from 



the ergosphere can be expected in models which assume that 
strong release of the accretion energy occurs in a corona sur- 
rounding the disc. Although details of the heat deposition 
(presumably involving buoyancy and reconnection of mag- 
netic fields) are uncertain, the spatial distribution of the 
heating is most likely to follow the dissipation distribution 
in the disc. 

Simulations of accretion driven by turbulent stresses 
generated by the magnetorotational instability (e.g. Haw- 
ley & Balbus 2002; De Villiers et al. 2003) support a model 
with an inner corona formed within a few gravitational radii. 
These simulations commonly indicate formation of a hot in- 
ner torus, of the size slightly exceeding r ms , at the inner 
edge of a nearly Keplerian disc. For a rapidly rotating black 
hole the inner torus, which can be plausibly related to coro- 
nal activity of the inner region, is located mainly within the 
ergospheric region. 

Observations of black hole systems indicate that at high 
accretion rates geometry of the innermost region is similar 
to that considered in this paper. High temperature of ther- 
mal component, dominating the soft-state spectra of black 
hole binaries, implies disc extending down close to a black 
hole (e.g. Gierlihski & Done 2004). A similar conclusion 
comes from the amount of X-rays reflected from the disc 
in a Seyfert galaxy, MCG-6-30-15 (Lee et al. 1999), which 
seems to be a high accretion rate object (Nowak & Chiang 
2000). Furthermore, Fe Ka line profiles indicate a very cen- 
tral location of the hard X-ray source in high accretion rate 
objects, including XTE J1650-500 (Miller et al. 2002) and 
GX 339-4 (Miller et al. 2004) in their soft states as well as 
MCG-6-30-15 (Wilms et al. 2001). In our subsequent paper 
we find that relativistic line profiles observed in black hole 
systems can be explained in a model involving illumination 
of a disc by a central corona and that the extreme profiles 
noted above require such a corona to be constrained within 
the ergosphere. 

Given the above observational evidence, we expect that 
results of this paper should be relevant to modeling observa- 
tions of objects with high luminosity to Eddington luminos- 
ity ratio. On the other hand, we emphasize that nonthermal 
tails observed in some soft state spectra (e.g. McConnell et 
al. 2002) cannot be due to special or general relativistic ef- 
fects in a thermal plasma. Such a tail most likely indicates 
that Comptonizing plasma contains nonthermal (in addi- 
tion to thermalized) electrons (e.g. Zdziarski 2000), which 
has not been taken into account in this study. 

Observations in lower luminosity states, including most 
normal Seyfert 1 galaxies and black hole binaries in their 
hard states, require that the innermost part of accretion 
disc is replaced by an optically thin flow, with the opti- 
cally thick material extending down only to 20 — 50_R g (e.g. 
Done, Madejski, Zycki 2000; Miller et al. 2001; Zycki, Done 
& Smith 1998; 1999). Hard X-ray spectra of such objects 
are well fitted by thermal Comptonization models, usually 
using a uniform temperature and density distribution, with 
electron temperature within 50-100 keV and a Thomson op- 
tical depth close to unity (e.g. Zdziarski 2000). If the spacial 
extent of such a uniform, hot plasma corresponds to the dis- 
tance of truncation of the disc, then the relativistic effects 
considered in this paper should not be important for the 
spectra. 

However, if the plasma is heated due to dissipation of 
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gravitational energy in a hot inner flow around a rotating 
black hole, then enhanced emission from the innermost re- 
gions can be expected. Accretion flows around rapidly rotat- 
ing black holes release a large fraction of the gravitational 
power within the ergosphere, as noted above. Furthermore, 
advection increases electron temperature in the innermost 
region of optically thin flows (Esin et al. 1997). Then, emis- 
sion of an inner corona may be very centrally concentrated 
and a detailed modeling is required to check importance of 
the gravitational effects for formation of spectra in such a 
case. 

Finally, we indicate that hints for orientation-dependent 
hard X-ray spectra, consistent with predictions of our model, 
come from recent studies of X-ray observations of Seyfert 
galaxies. In their analysis of OSSE observations, Zdziarski, 
Poutanen & Johnson (2000) found that Seyfert 2 galaxies 
have significantly harder spectra than Seyfert Is, with the 
average photon spectral index, in the 50 - 200 keV range, 
being approximately 2 and 2.5, respectively. As only in one 
Seyfert 2, considered in that paper, absorption of X-rays is 
Thomson thick, this result may indicate an intrinsic differ- 
ence between the hard X-ray spectra of Seyfert Is and 2s. 
A similar conclusion comes from BeppoSAX observations 
(Deluit & Courvoisier 2003), which indicate also that the 
cut-off energy is higher in Seyfert 2s than in Seyfert Is. 

According to the unification model of active galactic nu- 
clei (e.g. Antonucci 1993) the difference between type 1 and 
type 2 galaxies is only due to the viewing angle. Namely, the 
same object would be classified as a Seyfert 1 or as a Seyfert 
2 when observed at low and high inclinations, respectively. 
Following this unification scheme, we indicate that Seyfert 
2 galaxies may indeed be intrinsically harder than Seyfert 
Is due to gravitational focusing of X-rays by rotating black 
holes. 



6 SUMMARY 

We have developed the Monte Carlo model for studying 
Comptonization spectra emerging from accretion flows in 
the Kerr metric. We find that strong gravity effects are cru- 
cial for the formation of the spectra only in the ergosphere 
of a rotating black hole. Spectra involving a significant con- 
tribution from the ergosphere are strongly dependent on the 
observer's inclination. Furthermore, intrinsic curvature and 
spectral breaks, related to strong gravity effects, may occur. 

The high energy part of the ergospheric Comptonization 
spectrum is significantly attenuated when optically thick 
material is present in the equatorial plane of the Kerr ge- 
ometry. Therefore, a noticeable Comptonization radiation 
from the ergosphere is more likely in models with a tenuous 
plasma occupying the innermost region. 

We indicate that hardening of Comptonization spec- 
tra with increasing inclination, which is predicted by Comp- 
tonization in the Kerr metric, can be considered as an ex- 
planation for the difference between the hard X-ray spectra 
of Seyfert 1 and Seyfert 2 galaxies. 
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